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1. In curve fitting problems that arise in science or 
engineering we often demand that our approximating function be 
shape preserving in the sense that the approximation is convex 
when the data are convex or that the approximation is monotone 
and convex when the data are monotone and convex. In this paper 

we construct such shape preserving approximations by interpolating 
the data with polynomial splines of arbitrary degree. We 
formulate a regularity condition on the data which insures the 
existence of such a shape preserving spline, we present an 
algorithm for its construction, and we bound the uniform norm 
of the error which results when the algorithm is used to produce 
an approximation to a given f e C[a,b]. 

2. Let A :a=x A <x, <...<x = b denote an arbitrary 

n 0 1 n 

but fixed partition of the interval [a,b] with knots x^, let 
h. = X.,, - x., i = 0,l,...,n-l, and let h = max h.. Let 
Sp(k, £, A ) denote the space of polynomial splines of degree k 
and deficiency k - £ and assume that k = 3,4,... and £ = 1,..., 

I(k - 13/2] . 


*0n leave from the Institute of Computer Science, University of 
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Given real data (Xq,£q), (x ,f n ) we see ^ s £ Sp(k, l, A n ) 

such that 

(P) s(x^) = f^, i = 0,1,..., n and s is convex on [a,b]. 

The existence of a solution to (P) depends on the values of 
o i = (f i+1 - £ i )/h i , i = 0, 1 , . . . ,n-l . 

To construct s we set 

(1) s(x) = s i (x) , x i < x < x. +1 , i = 0 , 1 , . . . ,n- 1 
where 

(2) s i (x) = f.<j) 0 Ct) + £ i+1 <J) 1 (t) + h i [p i T 0 (t) - P i+1 T 1 (t)], 

t = (x - x i )/h i , 

where Pq»P]_ » • • • »P n are parameters to be determined later and 

1 t 

4> 0 (t) = X M k m (s)ds, * 1 (t) = J* q M k m (s)ds, 

Vt) = M k m 1 (s)ds, ^(t) = (1 - t)/ M k m ‘ 1 1 (s)ds, 

where denote the B-spline of degree i - 1 with knots 

0 = t n = ... = t - < t.., = ... = t. = 1. Clearly s e Sp(k, 5- , A ) 

U jj+i l J n 

with H = min{k-m-l, m} £ [ (k - l)/2]. 

Theorem 1 . The spline function s given by (1) - (2) solves 
the problem (P) if and only if 


(k - m)p i + m p i+1 


1 ^ i 


(k-m-l)p. + (m + l)p i+1 


for i = 0,1,..., n-1. If in addition p^ _> 0, then s is also 
nondecreasing on [a,b] . 


52 



3. The parameters which appear in (3) can be constructed 
as follows. We initially define 


l. = 

l 


ka . - (m + l)a . , . 
i v J i + . 


;-m- 1 


> i 0 , . . . , n- 2 , ^ ~ co y 


ka i 1 " O-m-l^ 1 

u o = -*» u i = irnr » i = 1 »--*» n » 

ka. - (k - m)p. 

v 0 = v i ' m ■ 1 ‘ 1 

and then for each i = 0,1,..., n set c. = max{u., £.}, 

’ ’ l i' l 

d. = min{v. , a . } . 
i i i 

Theorem 2. Let a. > a n and let a = +». If 
1—0 n 

(m + 1) O i+1 - a i ) > (k-m-l)(a i - a-_ 1 ) , i = l,...,n-2, 
then c^ £ d^, i = 0,...,n and any choice of p^ e [c^, d^] , 
i = 0,...,n makes (1) - (2) a solution of (P) . 

4. The error which results when the above procedure is 
used to construct a polynomial spline approximation to a continuous 
function f may be bounded as follows: 

Theorem 3 . Let f e C[a,b] be nondecreasing and convex on 
[a,b], and let s be the spline interpolant produced by using 


(1) - (3) with f^ = f(x^), i = 0,...,n and with 


11 ^0 - 


> 0 . Then 


U £ ‘ S L i Y^m^ 10 

where oo C f ; * 3 is the modulus of continuity of f and 

Yv m = 1 + A*k/m and i - ( k £ )2 k < A < 1 
k,m 2m' — — 


Moreover, for k = 2q - 1 and m = q - 1 (q > 1) 


4 - y 3,1 < y 5,2 < ‘ • - 2 ‘ 
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